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Abstract: We derive a formula for the light field of a monochromatic 
plane wave that is truncated and reflected by a spherical mirror Within 
the scalar field approximation, our formula is valid even for deep mirrors, 
where the aperture radius approaches the radius of curvature. We apply 
this result to micro-fabricated mirrors whose size scales are in the range of 
tens to hundreds of wavelengths, and show that sub-wavelength focusing 
(full-width at half-maximum intensity) can be achieved. This opens up the 
possibility of scalable arrays of tightly focused optical dipole traps without 
the need for high-performance optical systems. 
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1. Introduction 

Spherical mirrors are widely used, for example in the telescope, to collect light from a distant 
object and focus it to a point. In practice this produces not a point but a distribution of light over 
a finite region in the vicinity of the focus. Analysis of this behavior involves two classic topics in 
optics — diffraction from a circular aperture and spherical aberrations of the mirror. Typically, 
the spot size is characterized by a figure of merit such as transverse aberration or root-mean- 
squared blur radius. However, in a variety of modern applications, such as data storage, optical 
tweezers or atom trapping, it is only important to have a narrow central spot even though this 
may be accompanied by a broad distribution of low intensity in the wings. For such applications, 
the analysis requires a new approach, which we develop here. 

Our own interest is in making an array of atom traps so small that each trap can hold one 
atom (or none), but not two. If a second atom is caught, the two interact through a light-induced, 
inelastic collision and both are kicked out of the trap — an effect known as collisional blockade 
lUllll. The ejection of atom pairs occurs promptly provided the spot size of the trapping Ught 
is below ~ IjUm |3l, which is roughly the size of the optical wavelength A. One important ap- 
plication for such traps is the production of single photons on demand |l4l|5l|6l, which can then 
be used as a powerful resource for quantum cryptography and quantum information processing 

Qllllllol. 

There is a growing Uterature concerned with strong focusing of light onto single atoms in free 
space, as a means of achieving efficient atom-light coupling. A scheme proposed in Ref. IfTTI 
aims to emulate the field radiated from a single atom. This ideal can be approached using 
a parabolic reflector to focus light with highly sculpted intensity and phase profiles, but the 
scheme does not seem suitable for scaling to large arrays of traps. In Ref. lfT2ll . van Enk and 
Kimble considered the simpler case of focusing a uniformly polarized Gaussian beam with 
a lens that imposes a Gaussian phase profile. They were able to achieve a three-dimensional 
solution for the light field in the region of the focal point. Tey et al. recently extended this 
analysis to the case of a spherically converging Gaussian wave flSl , which gave somewhat 
better coupling. Both of these calculations neglect the aberrations produced by a real imaging 
system. 

In the work presented here, we consider focusing by spherical mirrors and we treat the aber- 
rations of the mirror exactly. No high-performance optical system is required and the method 
offers a simple way to achieve large arrays of traps that are tight enough to produce a strong 
collisional blockade for alkali atoms. The idea is to use arrays of hemispherical mirrors, typ- 




Fig. 1. (Color online) Light incident along the path r' is reflected at a point on the spherical 
mirror surface, then propagates along r to the point of observation (cartesian coordinates 



ically with radius of curvature R ~ 100 /im, that are lithographically etched directly into the 
surface of a silicon wafer lfT4l [T5]| and covered with a reflective coating. An incoming plane 
wave may then be focused to an array of points in a plane near the surface of the wafer, to 
produce a tight optical trap above each reflector. Once integrated into an atom chip |fT6l , these 
traps may be loaded either directly from a surface magneto-optical trap, or from a magnetic 
guide with optical molasses added to enable the blockade. 

In calculating the characteristics of these microtraps, we found it necessary to consider both 
the diffraction and the spherical aberrations. This led us to derive an analytical expression for 
the distribution of the light near the optic axis of a spherical mirror, which is applicable for 
wide-angle reflectors and does not seem to have appeared previously in the literature. Although 
the full problem can be solved numerically with a high-performance computer and sufficient 
patience, our formula has the advantage that it can be evaluated more or less instantaneously. In 
the rest of this paper, we derive our analytical results, compare them with full numerical simu- 
lations, and show that the silicon micro-mirrors are indeed suitable for making sub-wavelength 
atom traps. 

2. Theory 

Figure [T] shows a schematic micro-mirror with radius of curvature R centered on the point 
(0,0,0), illuminated by a point source at position (0,0, Zs)- We take zs so that the 

incident light is collimated and parallel to the optic axis. If the light were incident on a simple 
circular hole of radius a, the amplitude of the output in the far field propagating at angle a to 
the optic axis would be given by Ji [ka sin a)/ sin a, where k is the wavenumber 2k /X and Ji 
is the Bessel function ifTTl . As a starting point, let us consider ideal imaging that satisfies the 
Abbe sine criterion and focuses incident collimated light onto the plane z = R/2. In the image 
plane, the rays incident at angle a have an off-axis displacement of (/?/2)tana. The resulting 
Fraunhofer diffraction pattern in the focal plane has a radius at half-maximum intensity, r^^J^, 



given by 




diff 



(1) 



Here, we have neglected the difference between sin a and tan a, since we are only interested 
in apertures having a ^ A, which produce small angles of diffraction. The aperture size is 
conveniently expressed by the dimensionless ratio p = a/R. For a spherical mirror, p < 1, but 
for imaging that satisfies the Abbe criterion, the numerical aperture reaches its maximum value 
of 1 when the input radius is equal to the focal length (see ifTTl . §4.5.1). Thus the maximum 
aperture for our ideal optic is p = 1/2, giving a minimum spot size according to Eq. ([T]i of 
0.26 X. Throughout this paper we use the phrase "spot size" to mean radius at half-maximum 
intensity. 

This result captures the essence of the diffraction but does not treat the focusing correctly 
since a large-aperture mirror exhibits numerous orders of aberration. We therefore adopt Kirch- 
hoff's diffraction theory in order to estimate the spot size more accurately. The incident light 
is assumed to illuminate the mirror surface uniformly. Secondary waves then radiate to the ob- 
servation point xq = i^R, TjR, ^R), illustrated in Fig.[T] For simplicity we ignore the incoming 
beam and reflections from the flat part of the mirror because we are interested in the focal re- 
gion, where the reflected intensity is relatively high. Kirchhoff 's integral for the diffracted wave 
xj/f ^ is then 



— ^nCW-GWVn- 



^/jz/ , (2) 



where d^/ = R^d(j) dcos is an element of area on the curved surface of the mirror. The scalar 
r' is the modulus of the vector r' shown in Fig.[T] and similarly for r. The free-space Green 
function for the light field is G{r) = exp{ikr) / r, and Vn denotes the gradient taken along n, the 
normal to the mirror surface. Since the observation points of interest will be many wavelengths 
away from the mirror surface, we can assume that kr ^ I and therefore 

g'kr gikr 

Vn — « /^cos(n,r) — , (3) 
r r 

where cos(n, r) is the cosine of the angle between r and n. A similar expression holds for r — > r'. 

Taking zs ^ the source term exp(/A:r') /r' ^ \f/exp{ikRcos 0), with the complex ampli- 
tude y/ being constant. Then 



x^rf) = -J^jJ ; [cos(n,r)-cos(n,r')]^f^. (4) 



With this collimated input, cos(n,r') is just cos 6, while cos(n,r) is given by 

rcos(n,r) = R (1 — cos sin — i^cosQ) , (5) 

where we have taken Tj = without loss of generality. Since the light is to be tightly focused 
near the optical axis, we need only consider observation points where ^ <C 1 . Expanding r to 
first order around ^ = 0, we obtain 

n ?T — t: ^ sin 9 cos d 

= ^ + ^ 5i cos ^ . (6) 

Returning now to Eq. dUi, we eliminate the two cosines and replace rhy R5o everywhere except 
in the exponent. There, we keep theR^ 5i term as well and integrate over (j), giving 

(0) kxj/ /■/■e''<-(M+*+?5,cos,(.) / , . 

= ^ijj — ^ — 
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Fig. 2. (Color online) Intensity distribution on axis with R = lOOA for values of mirror 
aperture p = 0.2 (bottom), 0.4 (middle), 0.6 (top). These are offset vertically for clarity. 
Labels give the peak intensity for unit incident intensity, (a) Prediction of Eq. Q. Curves 
are re-scaled to be equal in height, (b) Result of full numerical integration of Maxwell's 
equations. Each curve has the same scale as the corresponding curve in (a). Dashed lines 
mark the geometrical focus at = 1 /2. 



where k = kR, jx — cos 9 and jXQ — {\— p^) This is our main result. 

In the hmit of small mirror aperture, i.e. small p and 0,6o^(l — C)'M^1'(1 — M^)'''^ 9 
and c/jU -Odd. Then 

Itf'hM(f) 

Thus the radius at half-maximum intensity is 

0/2 = 0.26(1 -C)A/p (whenp<l). (9) 

Evaluated in the focal plane ^ = 1/2, this duplicates Eq. ([Til, as one would expect in this limit 
where the mirror becomes an ideal thin element. The position of maximum intensity lies on the 
z axis, but not at = 1 /2. In this small aperture limit, the intensity increases as the mirror is 
approached because the spreading of the light due to diffraction dominates over the focusing 
due to the curvature of the mirror surface: this is the limit of small Fresnel number. 

3. Results 

Figure|2a) shows the intensity on axis versus position, as given by Eq. (|7]) for three larger values 
of the aperture. These all have Fresnel numbers greater than unity and exhibit peaks close to 
^ = 0.5. At first, as the aperture is increased, the light becomes more concentrated axially, 
giving rise to a narrower peak at p = 0.4 (middle) than at p = 0.2 (bottom). With even larger 
aperture, however, secondary structure appears on the large- side of the main peak. This is due 



q = 



(8) 



Fig. 3. (Color online) Simulated intensity distributions of light outside concave spherical 
mirrors with radius of curvature R = lOOA, etched in a plane substrate. The field of view 
is lOA (vertical) x40A (horizontal), centered on z = R/2, with the mirror to the right. 
The calculations are done by numerical integration of Maxwell's equations. Upper image: 
p = 0.1. Lower image: p = 0.4. 

to spherical aberration, i.e. to rays that are incident increasingly far from the axis and therefore 
cross the axis at larger values of ^ after reflection. In order to test these detailed predictions 
of Eq. (|7]i, we have integrated Maxwell's equations numerically using freely available software 
ifTSl . based on the finite-difference time-domain method |[T9l . with sub-pixel smoothing for 
accuracy at sharp interfaces EOl . Our resolution varies from 20 to 32 pixels per wavelength 
and we treat the mirror surface as a perfect conductor We exploit the cylindrical symmetry 
of the problem by taking the incident plane-wave beam to be circularly polarized. Figure [3] 
illustrates the solutions obtained in this way for p = 0.1 (upper frame) and p = 0.4 (lower 
frame). The numerical results for the intensity on axis are plotted in Fig. [2jb) for comparison 
with the plots shown in Fig. [2a). 

For the smallest aperture, there is good agreement on the shape and intensity of the curve, but 
the numerical integration exhibits additional rapidly oscillating fringes, similar to those seen in 
the upper frame of Fig. [3] These are due mainly to interference between the incident plane- 
wave and the field reflected from the mirror The same effect can be captured qualitatively by 
adding the incident field to Eq. (|7]i. Such an intensity distribution could be useful for making a 
tightly confining optical lattice but would not be chosen when a single, well-defined atom trap is 
required. The full numerical solution with medium aperture, shown in the lower frame of Fig. [3] 
and in the central curve of Fig.|2b), agrees very well with our analytical result. The interference 
fringes are much less evident here, as the amplitude of the reflected wave is much higher. There 
is also good agreement in the case of the largest aperture, except that the subsidiary structure 
due to spherical aberration is slightly smaller in the exact solution. This is because the electric 
fields of rays coming from the outer parts of the mirror are significantly inclined and should be 
added as vectors, whereas Eq. (|7]i, based on the scalar wave equation, adds them as scalars. We 
conclude that an aperture p = 0.4 is a good choice for achieving a single optical dipole trap 
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Fig. 4. (Color online) Normalized spot size r^j^j'^ as a function of aperture p with 
R = 100 A.. Solid line: radius given by Eq. Q evaluated at the peak of the axial inten- 
sity distribution. Dashed line: radius given by Eq. lUll for an ideal optic with focal length 
R/2. Dots: full numerical integration of Maxwell's equations. 



with tight axial confinement, and that this is well described by Eq. (|7]i. Note that the use of a 
larger aperture does not add substantially to the peak intensity. 

We turn now to the radial width of the intensity distribution. The dashed curve in Fig.|4]shows 
the ideal value r^^^, given by Eq. plotted against the aperture p up to the maximum value 
possible, p = 1/2. This suggests that apertures p > 0.2 should produce very small spot sizes. 
The solid line, derived from Eq. (|7]i, shows that although the spot size is somewhat increased by 
the mirror depth and the aberrations neglected in Eq. ([T]), it is still well below X /2. We calculate 
this width not in the focal plane ^ = 1/2 but at the value of ^ where the on-axis intensity is 
maximum, since that is where the optical dipole trap is actually formed. The dots show the 
widths obtained from the full numerical integration of Maxwell's equations. Near p = 0.2, 
these agree closely with the widths derived from Eq. (|7]i, but as the aperture opens, we see 
that the full solution gives a slightly larger spot size. This broadening is another manifestation 
of the vector nature of the light field. We conclude that parallel light, incident on a spherical 
micro-mirror with aperture p = 0.4, can produce an optical dipole trap that is well described 
by Eq. (|7]i, with spot size < A/2. 

Concerning the vector nature of the field, this is transverse on axis and follows the polarisa- 
tion of the input Ught. Off axis, however, the field acquires a z-component, given in the case of 
small p by 11111221 

,(0), 



(10) 



This produces the increase in ri/2 that we have noted above. We can estimate the size of the 
z-component from the derivative of Eq. (O, even when the mirror is not strictly in the small p 

limit. Near the focal plane and to first order in ^ , this gives xi/i^"^ / y/^f^ ~27rp^.j;/A. Forp =0.4 

this ratio grows from zero on axis to 0.4 at x = ri/2 = 0.4 A. Thus, the trapped atom will not 
see a constant polarization as it explores the volume of the trap, but assuming that its kinetic 
energy is significantly less than half the trap depth, the polarization will be approximately that 



of the input light. 




Fig. 5. (Color online) Normalized spot volume, as defined in the text, for varying aperture 
sizes and/? = 100 A. Dots: numerical integration of Maxwell's equations. Line: results from 
Eq. 0. 

For colHsional blockade experiments, the focal spot must have a small volume, in addition to 
small ri/2, in order to ensure no more than a single atom per trap. We define the spot volume 
Vi^2 the volume within the contour of half-maximum intensity, which is plotted in Figure |5] 
as a function of p. For small apertures, we find °^ '"1/2/'^' expected from dimensional 
arguments. For larger apertures, Vj /2 drops to a minimum near p =0.5, then rises to a level just 
below A^. As a point of reference, the spot volume for an ideal Gaussian TEMqo beam is given 
by ^1/2 = 41.8 r^^ij/X, which is 2.6 for ri/2 = A/2. Once again, Eq. (|7]i provides a very good 
approximation to the exact solution. Finally we note that although all the simulations shown in 
this article have taken R~ 100 A, this is not a critical requirement. When R is increased from 
60 A to 140 A, the spot volume obtained by numerical integration for p = 0.4 grows by less than 
14%, though of course the constraint on surface quality of the mirror becomes more demanding. 

4. Conclusion 

In conclusion we have derived a useful formula for the field distribution in front of a spherical 
mirror illuminated by a plane wave. Using this result we have shown that a modern micro- 
fabricated mirror of order 10-100 A in size can produce either a tight optical lattice, or a single 
spot with little additional structure, having a radius at half-maximum intensity well below A/2. 
The fabrication method makes it straightforward to scale this up to a large number of spots in 
any desired array. This method is therefore suitable for building atom traps for applications in 
quantum information processing, or for any other application where tight focusing is required 
and the low-intensity wings of the spot are unimportant. 
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